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ABSTRACT 

In order to better understand the nature of active region outflows, the electron 
density was measured by using a density-sensitive line pair Fe xiv 264.78A/274.20A. 
Since coronal line profiles of the outflow region are composed of a major component with 
a Doppler shift of < 10kms“^ and a minor component (enhanced blue wing; EBW) 
blueshifted by up to 100kms“^, we extracted EBW from the line profiles through 
double-Gaussian fitting. We tried applying the simultaneous fitting to those two Fe xiv 
lines with several physical restrictions. Electron density for both components (uMajor 
and nEBW) respectively) was calculated by referring to the theoretical intensity ratio as 
a function of electron density as per the CHIANTI database. We studied six locations in 
the outflow regions around NOAA AR10978. The average electron density was UMajor = 
]^q9.16±o. 16 gj^-3 nEBW = cm“^. The magnitude relationship between 

^Major and uebw was opposite in the eastern and western outflow regions. The column 
depth was also calculated for each component, which leads to the result that the outflows 
possess only a small fraction (~ 0.1) in the eastern region, while they dominate over 
the major component in the line profiles by a factor of five in the western region. 
When taking into account the extending coronal structures, the western region can be 
thought to represent the mass leakage. In contrast, we suggest a possibility that the 
eastern region actually contributes to the mass supply to coronal loops. 

Subject headings: Sun: corona — Sun: transition region — Sun: UV radiation 


1. Introduction 

Spectral coverage sensitive to the coronal temperature and unprecedented high signal-to-noise 
ratio of Hinode/EIS enabled us to reveal the existence of upflows at the edge of active regions 
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([Doschek et alJl2008l : iHarra et alJl2008l l. These upflows have been called “active region (AR) out¬ 
flows”, and are considered to be ejected from the bottom of the corona. It has previously been 
confirmed that these ou tflows persist for several days in the images taken by X-Ray Telescope 
(XRT) onboard Hinode (Sakao et alJ 2007). Some authors interpreted AR ou tflows as the source 
of the solar wind ( Harra et al. 2008 : Baker et al. 2009 : Brooks &: Warren 2011 1. 


Doschek et al.l ( 2003 ) analyzed emission line profiles of Fe xii 195.12A and revealed that the 


outflows are observed at the dark region outside an active region core. A preliminary result from 
EIS has shown that there is a clear boundary between closed hot loops in th e AR core 3 x 10® K) 
and extended cool loops 1 x 10® K) where the blueshift was observed ( Del Zannall2008l l . The 
upflows were seen in the low density and low radiance area. Meanwhile, redshift was observed in 
the AR core for all emission lines (Fe viii-xv). This apparent lack of signatures of any upflows 
at active region cores was explained as a st rong major com ponent closer to the rest wavelength 
in line profiles hinders the signal of upflows ( DoschekI 120121 '). but it has not been proved yet. The 
magnetic configuration of the ou tflow region has been modeled by ma gnetic field extrapolation from 


the photospheric magnetogram (IHarra et al 


2008 


Baker et al.ll2009l ). and it was revealed that AR 


outflows emanate f rom the footpoints of extremely long coronal loops in the edge of an active region 
(jHarra et al.l 120081 1. Close investigation revealed that AR outflows are located near the footpoints 
of quasi separatrix layers (QSLs), which form s the changes of the connectivity of the m agnetic helds 
from closed coronal loops into open regions ((Baker et al.ll2009l : (Del Zanna et al.ll201lh . 


The velocity of the outflow lies within the range of a few tens up to ~ 100 k m s ^. These veloc - 


ities were derived by subtractin g the fitted single-G aussian from raw line profiles (jHara et al.ll2008l ). 
and by double-Gaussian fitting (|Brvans et al.ll2010ll . By using extrapolated magnetic fields, the ac- 
tual velocity was d erived from the Doppler measurement and found to have a speed of 60-125 km s“^ 
(jHarra et al.M2008l l. The upflow velocity of AR outflows increase s with the formation temperature 


which emission lines Si vii-Fe xv represent (jWarren et al 
hotter emission line as 5 -20 km s~^ for Fe xii (formed at ^ 

(formed at ~ 3 x 10® K) ( Del Zannall2008l ). The appearance of the blueshifted regions often seems 


201 11 ). The blueshift becomes larger in 
1 X 10® K) and 10-30 km s“^ for Fe xv 


to trace loop-like structures. How ever, it is not completely understood whether the AR out flows 


are related to fan loop structures ( Warren et al.ll201ll : iTian et al.ll201ll : iMcIntosh et al.ll2012l ). 


AR outflows are observed as an enhanced blue wing (EBW) component in emission line profiles 
of Fe xii-xv. By fitting the line profiles by a single Gaussian, it was revealed that there i s a neg¬ 
ative correlation between blueshifts and line widths (jPoschek et al.l 120081 : iHara et al.ll2008l l. which 
indicates the existence of an unresolved component in the blue wing emitted from the upflowing 
plasma. This EBW does not exceed the major component by ~ 25% in terms of the intensity 


( Doschekll2012l ) 


Previous observations have revealed properties of the outflow from the edge of active regions 
such as (1) location: less bright region outside the active region core, (2) magnetic topology: 
boundary between open magnetic fields and closed loops, and (3) velocity: reaching up to u ~ 
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lOOkms”^ in the coronal temperature. Although a number of observations have revealed those 
physical properties, there remains one missing quantity: the electron density of the outflow itself. 
The density of an outflow region de rived by using the line ratio of Fe xii 186.88A/195.12A was 
~ 7 X 10® cm“® (IDoschek et al. 2008), which is sli ghtly lower than the typical value in active region 
(ue > 10® cm). Recently. iBrooks fc: WarrenI (2013) carried out differential emission measure (DEM) 
analysis at the outflow regions. It was revealed that the properties of DEM and also the chemical 
abundance are rather close to those of the active region, from which the authors concluded that 
the outflowing plasma originates in the active region loops. The interchange recon nection was 


Del Zanna et al.l 1201 il l. 


considered to be a candidat e for accelerating the plasma into the outer atmosphere (jBaker et al 

2OO9I: 


The electron density of the outflow itself should help us to better understan d the nature of the 


outfl ows. However, there have been few intensive attempts to do so until present (jPatsourakos et al 


2 OI 4 I I. One point of view is that those outflows are directly linked to the coronal heating in such a 
wav t hat the outflowing pla sma hlls the outer atmosphere and form the corona (jPe Pontieu et al 


2009 


McIntosh et al.ll2012l ). The impulsive heating in a coro nal loop induces an upflow from it s 


footpoint, which may account for what we see as the outflow (jPel Zannal l2008l : iHara et al.ll2008l b 
Outflows can be also c aused by the sudden change of the pressure environment in a coronal loop 
( Bradshaw et al. 2011 1. 


A theoretical estimation was recently proposed in terms of the ratio of the electron density 
betw een major comp onent (uMajor) and EBW component (ueew) in coronal emission line pro¬ 
files I Klimchuk 2 OI 2 I 1 . It was shown t hat if the tips of the chromospheric spicules supply the 


coronal plasma (IDe Pontieu et al.l I 2 OIIII . that ratio (here after denoted as nEBw/^Major) takes 


a value o f an order of 10-100, w hile tiny impulsive heating (i.e., nanoflare) creates the ratio 
of 0.4-1 (jPatsourakos et al.ll2014ll . Thus, it was suggested that the ratio UEBw/^^Major can be 


Patsourakos et al. ( 

2014) showed that this ratio peaks at order of unity, and suggested that type 

II spicules ( 

De Pontieu et al. 

2007 

) cannot be the primary source of the coronal plasma. 


In this study, we used the spectroscopic data obtained with EIS onboard Hinode in order to 
measure the electron density of the outflows. As a line pair suitable for our purpose, Ee xiv 264.78A 
and 274.20A were chosen because (1) those emission lines have a distinct enhanced blue wing at 
the outflow region which leads to better signal-to-noise ratio, (2) they consist of relatively clean 
emission lines and their line wings in the shorter wavelength side do not overlap with other emission 
lines, different from the cases for Ee xii 186.88A/195.12A and Ee xiii 202.04A/203.83A, and (3) the 
Ee XIV line pair is sensitive to the density range of rie = 10®^^^ cm“® as shown in EigurelH which is 


wider than other line pairs. The analyzed active region was as the same one as iPatsourakos et al 
(I 2 OI 4 I I. and one of our advantages is the spatial information (i.e., east/west edges), which was not 
focused in their study. 


The following parts of this paper are structured as follows. Section [2] describes EIS observation 
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Fe XIV 264A/274A 


Si VII 274A/275A 




Fig. 1.— Theor etical line ratio calculated by CHIANTI database version 7 (jPere et al.l 119971 : 


Landi et al.lpnid b (a) Fe xiv 264.78A/274.20A. (b) Si Vii 274.18A/275.35A. 


and wavelength calibration. Density of the outflows is derived in Section [3l and the results will be 
shown in Section 01 We propose a new technique for line profile analysis (A-ne) in Section [5l We 
discuss the nature of the observed outflows in Section [6l Section [3 will provide a summary of this 
paper. Two appendices describe some details in our analysis. 


Observation and calibration 


In this study, we analyzed a raster scan obtained with Hinode/EIS, which observed active 
region NOAA AR10978 (hereafter AR10978) at the center of the solar disk. The scan with narrow 
1" slit started on 2007 December 11 00:24:16UT and ended at 04:47:29UT. Field of view (FOV) was 
256" X 256" and exposure time was 60 s. The EIS data was processed through the standard software 
which detects the cosmic ray hits on the CCD pixels, subtracts the dark current bias, and corrects 
DN at warm pixels. The DN is converted into the unit of intensity: ergcm“^ s“^ sr“^ A . This 
quantity should be called spectral intensity in the literature. However, we use the term intensity for 
the sake of simplicity. One complicated point in the calibration is the thermal drift of the projected 
location on the CCD pixels due to the orbit al motion of Hinod e. We calibrated the absolute 
wavelength through the method developed bv iKamio et al.l (1201(11 ). Since the relative position of 
two emission lines Fe xiv 264.78A and 274. 20A is the most important f actor in th i s ana lysis, we 
carried out relative wavelength calibration whose details are described in iKitagawal (|2013l b 
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Fig. 2.— Line profiles of the active region AR10978. (a) Context image of AR10978 obtained on 
2007 December 11 00:24:16“04:47:29UT. Intensity of Fe xiv 264.78A is shown. Boxes numbered as 
1 (red) and 2 (blue) respectively indicate a footpoint of core loops and the outflow region, (b) Fe 
XIV 264.78A spectra, (c) Fe xiv 274.20A spectra. In each panel, line profiles at the footpoint of core 
loops {red histogram) and at the outflow region analyzed here {blue histogram) are shown in the 
upper half. The spectra were normalized by their integration. Residuals from the single-Gaussian 
fitting of each histogram are shown in the lower half. The Green histogram in panel (c) shows 
estimated spectrum of Si vii 274.18A. 


3. Density diagnostics of upflows 


One of our main achievements is density measurement of AR outflows. Previous observa¬ 
tions have revealed that the density of the outflow region measured by using a line pair Fe xii 
I 86 . 88 A/I 95 .I 2 A indicates 7 X 10^ cm ^ which is close to that of coronal holes rather than that of 
active regions (jPoschek et al.ll2008ll . However, density of the outflow itself, measured by separating 
its component from the major component in line profiles, has not been investigated so far. 

There are three reasons for the difficulties in the analysis of spectroscopic data obtained by 
Hmode/EIS. Firstly, the signals from an upflow are detected as an enhanced blue wing (hereafter, 
EBW) component in emission line profiles. Examples are shown in Figure [2j In each panel, line 
profiles at the footpoint of a core loop {red histogram) and at the outflow region analyzed here {blue 
histogram) are shown in the upper half. Residuals from single-Gaussian fitting of eac h histogram 


are shown in the lower half, which is quite useful in detecting weak signals in line wing (iHara et al 
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20081 ). There is a significant enhancement at the blue wing (< —100 km s both in Fe XIV 


264. 78A and 274. 20A as shown by blue histograms. Green histogram in panel (c) shows estimated 
spectrum of Si vii 274. 18A which was subtracted in the density diagnostics described later. The 
EBW component is weak in most cases as seen in spectra indicated by the blue histograms shown 
in Figure [2j In addition, EBW component is significantly dominated by the strong component 
closer to the rest wavelength, which makes the analysis of upflows quite uncertain. 


Secondly, the density measurement of the outflow itself needs the accurate determination of 
the rest wavelengths of emission lines from which we fit the two emission lines simultaneously and 
deduce the intensity. This is often laborious because we do not have the absolute measure of the 
wavelength corresponding to each observational spectral pixels. 

Thirdly, density measurement needs at least two emission lines from the same ion {e.g., Fe xiv 
as used in this paper). This means that the two emission lines should be fitted simultaneously using 
same parameters such as Doppler velocity and line width. No previous studies on the outflows from 
the edge of active region have dealt with such fitting. 

Our procedure of density diagnostics is as follows: (1) integration of neighboring multiple 
pixels in order to reduce the noise, (2) determination of the wavelength position corresponding 
to the same Doppler velocity, (3) removal of blending Si vii 274.18A from Fe xiv 274.20A using 
Si VII 275. 35 A as a reference, (4) simultaneous fitting of Fe xiv 264.78A and 274.20A, and (5) 
density inversion using a theoretical curve from CHIANTI as a function of the intensity ratio. In 
the following sections, each procedure will be described in detail. 


3.1. Integration of observational pixels 


The outflows from the edge of active regions are usually detected as an EBW in emissi on line 
profiles. Its intensity does not exceed ~ 25% of that of the major component (jDoschekll2ni2l ). This 
makes analysis difficult since the photon noise of the major component affects the emission from 
EBW. In addition, the region where the outflows can be seen is usually dark (i.e., small signal-to- 
noise ratio). In order to improve the signal-to-noise ratio, we integrated over multiple observational 
pixels in space using a square box with the size of 5" x 5". A larger integration box generally 
results in better signal-to-noise ratio. However, we chose that particular size of integration box so 
as not to lose the information of the outflow region. In the integration, the pixels with instrumental 
problems {i.e., hot or bad pixels) were excluded. 


3.2. De-blending of Si vii from Fe xiv 274. 20A 

Fe XIV 274.20 potentially has a contribution from Si vii 274.18A, which may become significant 
in the vicinity of an active region because Si vii emission often comes from the footpoint of cool loops 
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extending from the edge of the active region. We need to subtract this blend from Fe xiv 274.20A. 
In this study, the spectrum of Si Vii 274. IsA was calculated by using the observed line profile of Si 
VII 275. 35 A which is known to be clean (i.e., without any significant blend). The intensity ratio of 
Si VII 274 .I 8 A/ 27 5 . 35 A is at most 0.25 as calculated from CHIANTI version 7 ( Dere et al. 1997 : 


Landi et al.l 120121 '). The value has a dependence in the density range 10® cm ® < ng < 10 


10 


cm 


-3 


and it varies 0.06-0.27 (monotonically increasing) as shown in Figure [TJ First we remove the 
blending Si vii 274.IsA for the case rig = 10® cm“® (Si vii electron density), and after that we 
considered three cases of the ratio corresponding to the density of 10®, 10®, and 10^®cm~®. In 
order to make our analysis more robust, we excluded the location where the estimated intensity of 
Si VII 274.IsA exceeds 5% of the Fe xiv intensity. Using the theoretical ratio, the intensity of Si vii 
275. 35 A was converted into that of Si vii 274.IsA. The spectrum of Si vii 275.35A was then placed 
at Si VII 274.IsA taking into account the shift of Si vii 275.35A from the rest wavelength using 
the relative difference between wavelength of Si vii 274.IsA and 275.35A (i.e., I.ISOSA) given by 
CHIANTI database. Note that since there were no locations where Si vii 274.IsA dominates Fe 
XIV 274.20A in the data, we could not determine the relative wavelength position of the two Si vii 
lines, therefore we used the wavelength difference given by CHIANTI for the Si vii lines. The data 
points of the estimated Si vii 274. IsA in the wavelength direction were interpolated into the data 
points of Fe xiv 274.20A by cubic spline. Thus, we removed the blended Si vii 274.IsA from Fe 

XIV 274.20A. 


Concerning Fe xiv 264.78A, there are two possible blend lines: Fe xi 264.77A and Fe xvi 
265.OlA. As for Fe xvi 265.OlA, it is sufficiently far enough from Fe xiv 264.78A in non-flare situa¬ 
tions. Moreover, estimated peak intensity of Fe xvi 265.OlA was around 100ergcm“^ s“^ sr“^ A ^ 
in the observed outflow regior0, which is no greater than the background level of Fe xiv 264.78A 
as seen in Figure [3J Unfortunately, our data set did not have any isolated Fe xi emission line, 
which makes difficult to remove the blending Fe xi from Fe xiv 264.78A. Nevertheless, our crude 
estimation of the intensity of Fe xi 264.77A from Fe xi 188.21/188.30A (7254 77^/7i88 21 A ~ 0.03) 
leads to the potential influence on Fe xiv 264.78A by up to 5% in maximum. It is inferred from 
Appendix iBl that the error in our results can be considered to lie within ~ 3%. 


3.3. Simultaneous fitting of the two Fe xiv emission lines 

In order to make the htting more robust, the two emission line prohles of Fe xiv 264.78A/274.20A 
were fitted simultaneously. It is based on the consideration that the emission line profiles coming 
from the same ion species must have the same Doppler shift and the same Doppler width. As seen 
in Figure [21 emission line profiles of Fe xiv 264.78A and 274.20A from the active region core {red 


^We estimated the intensity from Fe xvi 262.98A included in EIS data. The line ratio Fe xvi 265.01A/262.98A 
was determined in the raster scan which started from 10:25:42UT since it included the spectra of both Fe xvi 262.98A 
and 265. OlA, and it resulted in the ratio of 0.083. 









(a) Equal width fitting (Fe XIV 264) 


(b) Equal width fitting (Fe XIV 274) 



264.6 264.7 264.8 264.9 274.0 274.1 274 . 2 _ 274.3 

Wavelength (A) Wavelength (A) 


(c) Independent fitting (Fe XIV 264) (d) Independent fitting (Fe XIV 274) 




(e) Simultaneous fitting (Fe XIV 264) 



(f) Simultaneous fitting (Fe XIV 274) 



Fig. 3.— Fitting results for Fe xiv 264.78A and 274. 20A obtained by three different models. Upper 
row (panels a and b): fitting with double Gaussians that have the same line width applied to each 
line profile independently (model 1). Middle row (panels c and d): fitting with double Gaussians 
that do not necessarily have the same line width applied to each line profile independently (model 2). 
Lower row (panels e and f): fitting with double Gaussians applied to two line profiles simultaneously 
without the assumption of the same line width of two components (model 3). 


histogram) are obviously symmetric, while those from the outflow region {blue histogram) have an 
EBW, from which it is not likly considered that a strong major component hinders any signals of 
the upflows in the active region core. This EBW did not exceed the major component anywhere 
in the outflow region (< 30%). Previous observatio ns have never s hown such emission line profiles 
whose EBW dominates over the major component (jDoschekll2012l ). 


In this study, the emission line profiles of the outflow region are assumed to be composed of 
two Gaussian components. Most previous analyses on the outflows at the edge of an active region 
assumed that the main component and EBW have the same line width in order to reduce avoid 
an unrealistic solution in the fitting parameter space, but the assumption could strongly affect the 
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fitting ([Bryans et alJ 120101 : iBrooks &: Warred 120121 ') . (Brooks &: WarrenI (|2012l i mentioned that this 
assumption may lead to the underestimation of the intensity of EBW. Line profile with EBW often 
shows rather longer tail in the line wing than could be represented by a Gaussian which has the 
same line width as the major component. Moreover, the assumption that the major component 
and EBW have the same line width is not based on the physical principles. 

In order to examine the differences in the fitting result between different constraint on the 
fitting parameters, we applied three fitting models to a line profile pair of Ee xiv 264.78A. Line 
centroid and line width are respectively denoted by A and W. The suffixes below represent: “1” 
for Ee XIV 264.78A, “2” for Fe xiv 274.20A followed by the component either “Major” or “EBW”. 
First model (model 1) assumes lEi^Major = kLi^EBW and W 2 , Major = kE2,EBW) and fits the line 
profiles of Fe xiv 264.78A and 274. 20A separately with double Gaussians that have the same line 
width for each component. The second model (model 2) also fits the line profiles of the two Fe xiv 
separately, but with double Gaussians that do not necessarily have the same line width for each 
component. The third model (model 3) fits the two Fe xiv line profiles simu ltaneously by applying 
A 2 , Major = aAi, Major, A2,EBW = aAl,EBW (« = 1.0355657, |Kitagawa||2013|) , VEi, Major = IT 2 , Major, 
and VEi^ebw = bL 2 ,EBW- We adopted model 3 for the electron density measurement in this study 
because it is physically most reasonable in the sense that the model calculates the parameters (line 
centroids and thermal widths) consistently for both emission lines and does not impose artificial 
restrictions on the line widths. 

The results for those three models are shown in Figure [3l We obtained smaller and more 
blueshifted s econd component (EBW ) with the model 1 in panels (a) and (b), which confirms the 
suggestion in IBrooks &: WarrenI (|2012l L In contrast, larger and less blueshifted EBWs were obtained 
with models 2 and 3 as clearly seen in panels (c)-(f). In addition to this, the line widths of EBW 
component were much broader for models 2 and 3 than for model 1. It is not clear at present 
whether the increased widths may indicate superposition of multiple upflow components, which 
will be another point to be revealed in the future. The comparison between those three models 
shows that the results in previous analyses probably underestimate the intensity of EBW with 
an artificial assumption that two components in line profile have the same line width. Moreover, 
independent fitting applied to two emission lines causes a discrepancy as seen in panel (c) and 
(d). The Doppler velocity of EBW component was —81.4kms“^ for Fe xiv 264.78A while it was 
—70.1kms“^ for Fe xiv 274. 20A. Note that the rest wavelengths were determined from a limb 
observation on 2007 December 6, so these Doppler velocities have an uncertainty of 10kms“^ at 
most. 


3.4. Density inversion 

Now the densities of EBW and the major component can be obtained by referring to the 
theoretical intensity ratio of Fe xiv 264.78A/274.20A as a function of electron density shown in 
Figure [TJ The intensity ratio monotonically increases within the density range of 10®cm“^ < Ue < 
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10^^ cm“^. The electron density in the solar corona generally falls between 10®cm“^ (for coronal 
holes) and 10 ^^ cm“^ (for flare loops), so the intensity ratio of Fe xiv 264.78A/274.20A is quite 
useful. The error in the density was calculated by using the l-cr error in the intensity ratio. The 
electron density is obtained from 


ne = F ^ 



( 1 ) 


where F~^ is the inverse function of the theoretical intensity ratio, and /264 and 1274 are respectively 
the observed intensity of Fe xiv 264. 78A and 274.20A. Using as the error of observed 

intensity ratio, we estimate the error of the density an^ as 


Ue ± iTn. = F~ 


-1 / -^264 

I274 


^ *^^ 264 /^ 74 ^ 


( 2 ) 


The error cr^e was not dealt symmetrically in this definition, which comes from the fact the function 
F has a curvature which can not be negligible compared to . 


4. Density derived from Fe xiv 264.78A/274.20A 
4.1. Results from single Gaussian fitting 


First we describe the results deduced from single-Gaussian fitting. As described above, line 
profiles at the outflow regions are known to have a distorted shape which cannot be well repre¬ 
sented by single Gaussian. Nonetheless, the results deduced from single-Gaussian fitting may be 
useful because the fitting is much more robust in terms of the freedom of variables {e.g, 4 pa¬ 
rameters for single Gaussian with constant background and 7 parameters for double Gaussians). 
Figure U] shows the map of intensity, Doppler velocity, line width of Fe xiv 264.78A, and elec¬ 
tron density derived from the line ratio Fe xiv 264.78A/274.20A. The blending Si vii 274.18A 
was taken into account and subtracted by referring to Si vii 275.SSA. It is clear from panel (b) 
that the outflow regions are present (i.e., blueshift) at the east/west edge of the active region core 
around {x,y) = (—280",—120") and (—175",—125"). Panel (c) shows that the line width at those 
outflow regions is larger than other locations by AW = 0.020-0.027A (square root of the differ¬ 
ence of squ ared line width) equ i valent to Sv = 20 -30 km s~^. which is similar to a result reported 
previously (iDoschek et al.l 120081 : iHara et al.l l2008l ) . The electron density at the outflow regions is 
He = 10®'^^®'® cm“^, which is lower than that at the core (rig > 10®'® cm“^). 


We defined the outflow regions as the locations (1) where the line width of Fe xiv 264. 78A is 
enhanced, and (2) which can be separated from fan loops seen in Si vii intensity map (not shown 
here). The selected six regions are indicated by white boxes in each map (numbered by U1-U6 
as written in panel (a), whose size is 8" x 8". Those regions are located beside the bright core as 
seen in the intensity map (panel a). We hereafter refer to U1-U2 as the eastern outflow region and 
U3-U6 as the western outflow region. 
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(a) Fe XIV 264 intensity 


(b) Fe XIV 264 Doppler velocity 
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(d) Fe XIV 264/274 electron density 
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Fig. 4. — Physical quantities deduced from single Gaussian fitting for Fe xiv 264.78A obtained on 
2007 December 11 00:24:16~04:47:29UT. (a) Intensity of Fe xiv 264.78A. (b) Doppler velocity of 
Fe XIV 264.78A. (c) Line width of Fe xiv 264.78A. (d) Electron density derived from the line ratio 
Fe XIV 264.78A/274.20A. 


4.2. Density of the upflows 

The electron density of EBW component was measured through the analysis described in Sec¬ 
tion O Eigure [S] shows the distributions of electron density for the major component (uMajor) in 
panel (a) and EBW component (nEsw) in panel (b). Pixels where the peak intensity of the major 
component (/Major) did not exceed 2.0 x 10^ ergcm“^ s“^ sr“^ A were masked by black. This 
threshold was determined by using the scatter plot of intensity and electron density of the ma- 
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(a) rig (major component) 



(b) Og (EBW component) 
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Fig. 5.— Electron density map deduced from two Gaussian fitting of an emission line pair Fe 
XIV 264.78A/274.20A obtained by the raster scan on 2007 December 11 00;24:16“04:47:29UT. (a) 
Electron density of the major component, (b) Electron density of EBW component. Same color 
contour are used in the two panels. Pixels where the peak intensity of the major component (/Major) 
did not exceed 2.0 x 10 ^ erg cm“^ s“^ sr“^ A were masked by black. White boxes numbered U1~U6 
are the same as those in Figure 01 The white dashed box indicate the entire western outflow region. 


jor component shown in Appendix O Pixels falling into the next three conditions were displayed, 
and others were masked by black. (1) /Major A 2.0 x 10^ ergcm“^ s“^ sr“^ A . (2) The intensity 
of EBW component (/ebw) exceeds 3% of that of the major component (/EBw//Major A 0.03). 
(3) The difference between the Doppler velocity of EBW component (uebw) and that of the ma¬ 
jor component (uMajor) satishes uebw “ ^Major < —30kms“^ {i.e., the two components are well 
separated). 


The relationship of electron density between the major component and EBW component are 
shown in Figure [6l Scatter plot in panel (a) shows the electron density for the outflow regions 
U1-U6 {colored symbols) and for the entire western outflow region indicated by the white dashed 
box in Figure [5] (WacA; dots). The eastern outflow regions (U1-U2) and west ones (U3-U6) exhibit 
different characteristics. The scatter plots for U1-U2 indicate nMajor < ^i-ebW) while those for U3- 
U6 indicate u-Major A it-ebw- Panels (b) and (c) show the same data but in histograms for which 
colors again indicate the selected outflow regions. The gray (the major component) and turquoise 
(EBW component) histograms in the background of panel (c) are made for the entire western 
outflow region. Those two histograms clearly indicate that uebw (10®'®^^*^'^^ cm“^) is smaller than 
ii-Major (10®'^®^*^'^^ cm“^) at the entire western outflow region, which confirms that our selection of 
the studied regions was not arbit rary. Note that our resul ts for the western outflow region (U3-U6) 
roughly consistent with those of iPatsourakos et al.l (j2014l i (i.e., nEBw/^Major ^ !)■ 
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Fig. 6.— (a) Scatter plot for Fe xiv electron density of the major component vs. that of EBW 
component. Colors show the selected region indicated by white boxes in Figure [H Triangles 
(Diamonds) represent the data points in the eastern (western) outflow regions. Numbers beside 
data points correspond to the name of the white boxes. Black dots show the electron density for the 
western outflow region indicated by a white dashed box in Figure [5l The dashed line indicates the 
point where two densities equal each other, (b) Histograms for the electron density of the major 
component (dotted) and EBW component (solid) in the eastern outflow region, (c) Histograms for 
the electron density of the major component (dotted) and EBW component (solid) in the western 
outflow region. The gray (the major component) and turquoise (EBW component) histograms in 
the background are made for the entire western outflow region. Those two histograms are multiplied 
by 0.1. 


4.3. Column depth 


Using the obtained electron density for each component in Ee xiv line profiles, the column 
depth of each component can be calculated. We use the equation for the column depth including 
the filling factor, 


h* = hf = -j- -- 

niG(ne, T) 


( 3 ) 


where / is the filling factor, I is the intensity of an emission line, Ug is the electron density, and 
G(ne, T) is the contribution function of an emission line. The quantity h* physically represents the 
plasma volume per unit area along the line of sight. Here the temperature substituted to Equation 
([3]) was simply assumed to take a single value Tf at which the contribution function G(ne, T) 
becomes maximum (logTf [K] = 6.30 for the Fe xiv lines used here). Panel (a) in Figure [ 7 ] shows 
a scatter plot for the column depth of the major component (/iMajor) and that of EBW component 
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Fig. 7. — (a) Scatter plot for column depth of the major component vs. that of EBW component. 
Colors show the selected region indicated by white boxes in Figure 01 Triangles [Diamonds) 
represent the data points in the eastern (western) outflow regions. Numbers beside data points 
correspond to the names of the white boxes. Black dots show the column depth for the entire 
western region indicated by a white dashed box in Figure O The dashed line indicates the point 
where two quantities equal each other, (b) Histograms for column depth of the major component 
[dotted) and FEW component [solid) in the eastern outflow region, (c) Histograms for column 
depth of the major component [dotted) and EBW component [solid) in the western outflow region. 
The gray (the major component) and turquoise (EBW component) histograms in the background 
are made for the entire western outflow region. Those two histograms are multiplied by 0.1. 

(^EBw)- Colored symbols respectively indicate the studied regions (U1~U2 for the eastern outflow 
region, and U3-U6 for the western outflow region). Similar to the result for the electron density, 
the eastern and western outflow regions exhibit different characteristics: /iMajor > ^ebw hr the 
eastern region, and /iMajor < ^ebw in the western region. Panels (b) and (c) display the same data 
in the form of histograms for the eastern and western outflow region respectively. The gray and 
turquoise histograms in the background of panel (c) show the results for the entire western outflow 
region indicated by a white dashed box in Figure [5j Table [T] shows the column depths averaged in 
each studied region. 

The result /iMajor < ^ebw in the western outflow regions (U3-U6) means that the upflow 
dominates over the major component in terms of the volume, opposite to the composition ratio 
of emission line profile itself. The value of /iebw — cm can be understood by considering 

that the inclination of the magnetic field lines in the western outflow region was 30°-50° (given 
the potential field calculation) and the horizontal spatial scale of the region was the order of 
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Table 1: Doppler velocities, electron densities, and column depths of EBW component and the major 
component derived through the double-Gaussian fitting applied to Fe xiv 264A/274A. Note that the 
Doppler velocities listed in the table are calculated by using limb spectra observed independently 
on 2007 December 6 as a reference of zero velocity, which leads to errors up to 10kms“^ at most 
mainly originating in the absolute wavelength calibration. 



EBW component 


The major component 

UDop(kms ^) 

log rie [cm 

log h [cm] 

UDop(kms ^) 

log Ue [cm 

log h [cm] 

Eastern outflow region 
U1 -92.4 ± 2.4 

U2 -84.8 ±21.4 

9.17 ±0.09 

8.95 ± 0.09 

7.03 ± 0.22 

7.67 ± 0.34 

-4.7 ±0.9 

-3.6 ± 1.7 

9.10 ±0.04 

8.93 ± 0.09 

8.38 ±0.05 

8.64 ±0.12 

Ave. 

-88.8 ± 15.2 

9.06 ±0.14 

7.36 ± 0.43 

-4.2 ±1.4 

9.01 ±0.11 

8.51 ±0.16 

Western outflow region 






U3 

-61.4 ± 15.7 

8.79 ±0.21 

8.25 ±0.35 

-0.6 ±2.6 

9.31 ±0.09 

7.74 ±0.12 

U4 

-56.3 ± 15.2 

8.53 ±0.17 

8.80 ± 0.34 

3.4 ±2.5 

9.34 ±0.10 

7.74 ±0.10 

U5 

-73.2 ± 12.2 

8.48 ±0.21 

8.59 ±0.41 

-1.3 ± 1.0 

9.11 ±0.07 

7.67 ±0.18 

U6 

-54.3 ± 13.3 

8.64 ±0.12 

8.40 ± 0.29 

-0.8 ± 1.1 

9.10 ±0.05 

7.89 ±0.10 

Ave. 

-62.0 ± 16.0 

8.60 ±0.22 

8.53 ±0.41 

0.1 ±2.7 

9.22 ±0.14 

7.74 ±0.15 


10'^(~ 10® cm), which leads to the vertical height of nearly the same amount. On the other hand, it 
is clearly indicated that /iebw is smaller than /iMajor by up to one order of magnitude in the eastern 
outflow region (U1--U2). This means that the upflows possess only a small fraction compared to 
the plasma characterized by the major component in line profiles. The Doppler velocities, derived 
electron densities, and the column depths for the studied outflow regions are listed in Table [H 

Note that in the line profile analysis, we assumed that the electron density corresponding to 
the temperature of Si vii (i.e., the transition region; hereafter nsivu) was 10®cm~^. We discuss 
this assumption and its influence on our results in Appendix iBl 


5. A-Ue diagram 

We modeled the spectra by the composition of two Gaussians in the above analysis. However, 
it is difficult to prove whether or not this assumption is suitable for the outflow regions. There are 
two alternative approaches to dealing with such a spectrum consisting of more than two Gaussians. 
One way is to adopt multiple-Gaussian functions (more than two components) and resolve multiple 
flows existing in a emission line. The more free parameters we use, the spectra would be fitted with 
less But this does not mean that we extracted a great deal of useful physical information from 
the spectra. The number of local minima increase with the complexity of the fitting model, and 
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the fitting process becomes an ill-posed problem. 


The other way is our new type of plot without assuming any htting model. Each spectral 
bin in a spectrum pair is used to derive electron density at each bin, which we refer to as “A-rie 
diagram” hereafter. With this method, we measured the electron density of the plasma which have 
the speed of UDop = c (A —Ao)/Ao (Aq: rest wavelength), which is a function of wavelength. Consider 
a density-sensitive pair of spectra (/>i(A) and 02(A) emitted from the same degree of an ion. These 
emission lines must have the same Doppler velocity because they came from the same degree of the 
ion, so after converting the variable A into Doppler velocity UDop as denoted by 0*(uDop) = 0i(A) 
(f = 1, 2), we can calculate the electron density as a function of the Doppler velocity 


nlivBop) = R ^ 


02('^Dop) 

_0l(^Dop) 


( 4 ) 


The derived n*(uDop) can be converted into a function of wavelength in either spectrum, ne(A), 
by the Doppler effect equation. Function R{ne) is the ratio of intensities from two emission lines 
which is a function of electron density, so when we know the intensities of two emission lines which 
are represented as 


h = J 0i(A)(iA, 

l2 = 1 02(A)dA, 

electron densities can be usually derived by 



( 5 ) 

( 6 ) 

( 7 ) 


Note that we used the same curve as shown in panel (a) in Fig. [T] for the function R{ne). This 
assumes that R{ne) is the same for all wavelengths in the range of interest, which we have not 
investigated in detail. 


As shown in the above equations, A-Ue diagram represents that of the particles which move 
with that speed, in other words, we do not obtain the electron density of the whole plasma as 
an ensemble of Maxwellian distribution. We emphasize that the advantage of our method using 
Equation Q is that even if we do not know the precise functional form of spectra, it gives us the 
electron density as a function of Doppler velocity without any modeling. 


5.1. Method 

Making A-Ug diagram contains the following processes: (1) subtraction of blending emission 
line, (2) adjusting wavelength scale of Fe XIV 264.78A to 274. 20A by interpolation, and (3) density 
inversion at each spectral pixel. Since the blend of an emission line Si Vii 274.IsA into Fe xiv 
274 . 20 A was already described in Section [321 here we explain only processes (2) and (3). 
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Since the EIS instrument does not have absolute wavelength scale, the corresponding wave¬ 
length location of the sam e velocity i n Fe x iv 264.78A and 274. 20A must be determined from the 
data itself as described in iKitaeawal ([2013). Using obtained relation Aobs,274/Aobs,264 = 1.03 5 5 6 57 
(±0.0000044), each wavelength value imposed on the spectral window of Fe xiv 264.78A was pro¬ 
jected onto the values on the spectral window of Fe xiv 274. 20A by the scaling 


A = aA264,i (a = 1.0355657), (8) 

where a number i indicates the ith spectral pixel in a spectrum of 264.78A. 

Since the wavelength values of each bin of projected Fe xiv 264.78 A do not generally coincide 
with those of Fe xiv 274. 20A, the projected spectrum was interpolated by a cubic spline in order 
to align two Fe xiv spectra in identical wavelength bins. 

We can calculate the ratio of spectral intensity Fe xiv 264.78A/274.20A at each spectral bin. 
Now we are able to derive the electron density in the same way described in Section 13.41 Because 
intensity at each spectral bin has larger errors compared to the integrated intensity {e.g., double- 
Gaussian fitting), the estimated errors for the electron density in the A-Ue diagram become large 
especially for the line wing. 


5.2. Verification of the method 


In order to test the validity of X-n^ method, we synthesized spectra of Fe xiv 264. 78A and 
274. 20A taking into account the spectral resolution of EIS and instrumental broadening. The 
spectra were composed of two components which represent plasma at the rest and an upflow. 
While the physical parameters for the major rest component (peak, Doppler velocity, and width) 
were hxed, those for a minor blueshifted component (i.e., upflow) were taken as variables. We 
made X-n^ diagrams for the minor component with 

• electron density of 8.50, 8.75, 9.00, 9.25, and 9.50 in the unit of log cm“^, 

• intensity of 1, 5, 10, 15, and 20% (ratio to the major component in Fe xiv 274. 20A), 

• Doppler velocity of 0, —50, —100, —150, and —200kms“^, 

• thermal width of 2.0, 2.5, 3.0, 3.5, and 4.0 MK. 


The nonthermal width was not considered in this test because essentially it does not produce any 
differences. In this paper, the tests only for electron density and Doppler velocity will be given 
below, since the dependence on them are significant. The other two variables (i.e., int ensity and 
therm al width) do not have strong effects and are described in the author’s PhD thesis ([Kitagawa 


2013). 
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5.2.1. Dependence on electron density 

The most important point on X-Ug diagram is whether it reflects the electron density of the 
components which compose spectrum properly or not. In order to test that, we synthesized the 
spectra which are composed of the major component at the rest which has the fixed electron 
density of logne[cm“^] = 9.0 and the minor component which has a variable electron density. Five 
cases (logUe [cm“^] = 8.50, 8.75, 9.00, 9.25, and 9.50) were analyzed, where the peak ratio of the 
minor/major component was 15% with fixed upflow speed v = —100kms“^. Panels (a) and (b) 
of Figure [8] show the spectra of Fe xiv 264.78A and 274. 20A respectively. Colors {blue, turquoise, 
yellow, green, and red) indicate the five cases calculated here. After converting the wavelength scale 
of 264.78A to 274. 20A, X-n^ were obtained as shown in panel (c) of Figure[8l The triangles in panel 
(c) indicate centroid and electron density of the given minor component. It is clear that those X-n^ 
diagrams clearly reflect the changes of the electron density from logUe [cm“^] = 8.50-9.50. Despite 
the spectra being composed of only two components, X-n^ diagrams do not become a step function 
but a smooth function. This is natural because the two Gaussians in the spectra contribute to 
each other by their overlapping wings. We claim that the method proposed here (A-Ue) is a good 
indicator of the electron density of components in the spectrum. 


5.2.2. Dependence on velocity 

The dependence of X-n^ diagram on the Doppler velocity of the minor component is obvious. 
The spectra of Fe xiv 264. 78A and 274. 20A, and X-n^ diagrams are shown in panels (d), (e), and (f) 
respectively of Figure[8l Colors indicate the hve cases for variable Doppler velocity calculated {blue: 

, turquoise: — 50kms“^, ^reen: —100kms“^, yellow: —150kms“^, and red: —200kms“^). 
Major rest component was at rest (0kms“^) with the electron density of logne[cm“^] = 9.0. 
The triangles in panel (f) indicate centroid and electron density of the given minor component. 
The relative intensity of the minor component is 15 % of that of the major component and the 
electron density of the minor component was set to log rie [cm“^] = 8.5 in all hve cases here. The 
location of dips in X-n^ diagram well represent the centroid position of the input minor component 
when two components are separated so that the spectrum is dominated by themselves near their 
centroids. This is not the case for v = —50kms“^ {i.e., yellow), where those two components 
are not separated so clearly. In this case, X-Ue diagram gradually decreases from longer to shorter 
wavelength. One advantage of the method described here is that we are able to know the tendency 
of electron density of uphow/downhow without any htting to the spectrum which might produce 
spurious results occasionally. 

The tests for the four variables {i.e., density, intensity, velocity, and thermal width) indicate 
that the method proposed here {X-Ug diagram) is a powerful diagnostic tool for coronal plasma which 
may be constituted of several component along the line of sight and form non-single-Gaussian line 
profile. In the next section, we exploit this X-Ug diagram so that the result obtained by double- 
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Fig. 8.— Upper panels: (a) synthetic line profiles of Fe xiv 264.78A, (b) those of Fe xiv 274. 20A, 
and (c) X-Ue diagrams. Each color indicates different electron density of the minor blueshifted 
component {blue: 8.50, turquoise: 8.75, green: 9.00, yellow: 9.25, and red: 9.25 in the unit of 
log cm“^). Electron density of the major component was fixed to logne[cm“^] = 9.00. The 
triangles in panel (c) indicate centroid and electron density of the given minor component. Lower 
panels: (d) synthetic line profiles of Fe xiv 264.78A, (e) those of Fe xiv 274. 20A, and (f) A- 
Ue diagrams. Each color indicates different velocity of the minor blueshifted component {blue: 
0km.s~^, turquoise: —50km.s~^, green: —100kms~^, yellow: —150kms“^, and red: —200kms“^). 
The major rest component was at rest (0kms“^). 


Gaussian fitting would be confirmed {i.e., upflows are more tenuous than the rest component). 
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Fig. 9.— Intensity map of Fe xiv 264. 78A obtained with EIS. Five arrays of colored diamonds 
[red-violet) are the locations where X-Ug diagrams were made. The locations cut across the active 
region core and the outflow region. Orange contour indicates the line width of 0.035A. 

5.3. X-Ue diagram in AR10978 

The electron density of the outflow region in AR 10978 is investigated through X-Ug diagram 
here. Figure [9] shows intensity map of Fe xiv 264. 78A obtained with EIS. Orange contours indicate 
the line width of 0.035A, which becomes an indication of the outflows. Five horizontal arrays of 
colored diamonds [red-violet) which cut across the active region core and the outflow region are 
the locations where we made X-Ug diagrams. First, we look at the location indicated by blaek plus 
signs named C (core) and U (outflow). 

In Figure [TUI the line profiles of Fe xiv 274. 20A, interpolated 264. 78A and estimated Si vii 
274. IsA (see Section 13.21) are shown by solid, dashed, and dotted spectrum respectively in panel (a) 
for the active region core and (b) for the outflow region. We can see an enhanced blue wing in line 
profiles of Fe xiv in the outflow region. The vertical dashed lines indicate rough reference of the 









- 21 - 



(c) A,-ng diagram at C (d) X-rig diagram at U 



Fig. 10.— (a) Line profiles of Fe xiv 274.20A (solid)/264.78A (dashed; adjusted to the wavelength 
scale of 274. 20A) and Si vii 274.IsA (dotted) estimated from 275.35A at the active region core, (b) 
Line profiles at the outflow region, (c) X-Ue diagram at the active region core, (d) A-rie diagram at 
the outflow region. 

rest wavelength position A = 274.195A which was the average line centroid above the limb in the 
2007 December 18 data (possible error up to O.OlA). 

Panels (c) and (d) in Figure [lOl show the X-Ue diagram for the active region core and the 
western outflow region respectively. The horizontal green dotted line in each plot indicates the 
electron density averaged in the neighboring three spectral bins which are nearest to A = 274. 20A 
(i.e., rest wavelength). Those X-n^ diagrams in the two locations exhibit a different behavior at the 
shorter wavelength side around A = 274.00-274. 20A: the diagram in the active region core is roughly 
constant while that in the western outflow region slightly decreases at the shorter wavelength. The 
number written in the upper left corner of each plot indicates the linear slope fitted within the 
wavelength range A < 274. 20A. This implies that the electron density of the outflows (i.e., shorter 
wavelength side) is smaller than that of the major component closer to the rest wavelength. 

In order to confirm the above implication more robustly, we see the variation of X-n^ diagram 
along X direction from the active region core to the outflow regions. The selected region spans from 
the active region core (red diamond) to the outflow region (violet diamond) as seen in Figure [9j 
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The boundary of the active region core corresponds to the color between yellow and light green. 
The A-Ue diagrams at each cut (1-5) are plotted in Figure [TTJ We can see clear changes of the 
A-rie diagrams with colors. The A-ng diagrams for cut 1 show a small hump around 274.00-274.lOA 
representing that EBW component has larger electron density than the major component, though 
the hump at almost all locations {red-blaek) might mean that it was caused by an anomalous pixel 
(e.^., warm pixel). Both for cut 1 and 2, the diagrams show flat or slightly decreasing behavior as 
a function of wavelength at all locations. These behaviors are consistent with the result obtained 
in Section [4] (region U1 and U2) which indicated that the electron density of the outflows in the 
eastern edge is almost the same or slightly larger. On the other hand, in the western outflow region 
(cut 3-5), those for the outflow region show a dip around 274.lOA. This wavelength corresponds to 
V = —llOkms”^ for the emission line Fe xiv 274.20A, from which it is implied that the outflows 
in the western edge are composed of less dense plasma compared to the plasma characterized by 
around 274. 20A existing along the line of sight. Note that this velocity does not mean that of the 
upflows because no fitting was applied in A-rig diagram. 

The electron density of EBW component evaluated from A-rig diagrams around A = 274. lOA 
was log rig [cm“^] = 9.0-9.2 in the eastern outflow region, and logUg [cm~^] = 8.5-9.0 in the western 
outflow region, which also coincides with the result obtained through the double-Gaussian fitting. 
By exploiting A-rig diagram as a new diagnostic tool, we can now support the results obtained in 
Section HI 


6. Discussion 


6.1. Theoretical estimation of electron density 


De Pontieu et al.l (|201lh proposed that the tip of the spicule is heated up to the coronal tem¬ 


perature (though the heating mechanism has not been revealed), and is injected to the higher 
atmosphere where the heated plasma form the coron a. The ele c tron d ensity of upflows from the 
tips of the spicules is estimated by Equation (10) in iKlimchukI (120121 1 which considers the mass 
conservation, 

nvp^sShs = UchcA, (9) 

where nup,s is the electron density of an upflow (a suffix s denotes spicule), 6 is the fraction of 
the spicule that is heated to coronal temperatures, hg is the height of the spicule, Uc is the coronal 
density after the tip of the spicule expands into the corona, is the length of coronal loops, and A 
is the expansion factor of the cross section of coronal loops from the chro mosphere to the corona . 
Using typical coronal values: Uc — 10®cm“^, he — 5 x 10® cm, (5 ~ 10% ( De Pontieu et al.ll201ll ). 
hg — 10® cm in the maximum height, and ^ ~ 10 (this factor has not been determined precisely 
yet, but is larger than unity), the electron density of upflows is estimated as 


n„P,, == 5 X 10“ 


hr. 


5 X 10® cm 


_ 6 _ 

or 


-1 


hg 


10 ® 


cm 


-1 


cm 


-3 


( 10 ) 














log Hg [cm'^] log [cm'^] log [cm' 
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Fig. 11.— X-Ue diagrams at the locations indicated by colored diamonds in Figured (Cuts 1 and 
2 ; including the eastern outflow region, Cuts 3-5; including the western outflow region). 

For impulsive heating, giving the typical energy content of nanoflares (ie., 10^^ erg) and con- 
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sidering the enthalpy flux as a response of the transition region below the corona leads to 


5 


Ei 


TTr-^Ti 


( 11 ) 


where p is the gas pressure of the upflow, uup,i is the speed of the upflow, Ei is the released 
energy by the impulsive heating, rgt is the radius of the coronal strand (ie., thin coronal loop 
as an elemental structure), and ti is the duration of the impulsive heating. Kinetic energy flux 
can be neglected because the upflow speed is around half the speed of sound (~ 200kms“^ at 
logT [K] = 6.3), which means the ratio of the kinetic energy flux to the enthalpy flux is the order 
of 0.1. Typical parameters Ei ~ 10^^ erg, uup,i — 100kms“^, rst ~ 10'^^®cm and ri ~ 10-100 s (this 
value contains a large degree of uncertainty because of a lack of knowledge at present) imply 


nup,i ~ 5 X 10 ' 


8-10 


Ei 


1024 gj.g J \ ]^q7 8 


rst 


-2 


10 s 
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-1 


t’UP,! 
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cms 


-1 


-1 


cm 


-3 


( 12 ) 


for which we used p = 2nup, i/cpTi wher e uttp, i is the electron d ensity of the upflow and Tj is its tem¬ 
perature. Recent observat ion bv Hi-C f Kobayashi et al.ll2014l ) indicated that the width of coronal 
strands is around 450km ( Brooks et al.ll2013l ). Note that Ei and Ti have not been observationally 
constrained well so far, though we used values which are considered to be reasonable at present. It is 
clear that the predicted electron density estimated by adopting the typical coronal values from the 
spicule and impulsive heating significantly exceed the derived upflow density (npEW < 10® cm“^ in 
our analysis). Equation (fT^ can be used to estimate the parameter range where the predicted up¬ 
flow density becomes similar to the observed value since there is much uncertainty in the parameter 
Ti. For example. If the heating continues for Ti = 500 s, Equation (I12jl leads to nup, i — lO^^^cm”^ 
(i.e., near the obtained upflow density) for the case other parameters keep their typical value. 


6.2. Mass transport by the outflow 


We estimate the mass flux of the outflowing plasma Tout in the western outflow region by 
using Doppler velocity and electron density of EBW component. The electron density was Ug — 

108.7 

cm and the Doppler velocity was —60kms 7, The total area (S) of the entire western 
outflow region was roughly 30" x 40" (S' ~ 6 x lO^^cm^). Considering the inclination angle of the 
magnetic held of 30°-50° as calculated by potential held extrapolation of an MDI magnetogram 
(mentioned in Section [6.3p . the speed of the outflow is roughly thought to be u ~ 70-90km s”^. 
Thus, Tout can be estimated as Tout = 2ngpvS = (4-5) x lO^® gs”^ where /r is a mean mass of ions 
which was set to 1 x 10“^^ g. For a comparison, we also evaluate the total mass contained in the 
active region. Using volume of U = (100")^ = 4 x 10^® cm^ and typical density rig = 10®^70gj^-3^ 
the total mass Mar is evaluated as Mar = 2ngpV = 8 x lO^'^^^s 


_ ^This implies that if the mass in the active region is actually lost by the outflow ([Brooks &: Warren 

2012l h the time scale of the mass drain becomes Tout = MAR/Eout = 2 x 10 ^^^ g (f.g.^ several hours 


to a couple of days). Since the lifetime of active regions is much longer than this time scale, up 
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Fig. 12.— EIS Fe xiv 264.78A intensity map. Orange (Turquoise) contours indicate a magnetic 
field strength of +250 (—250) G in an MDI magnetogram taken during the EIS scan. Six white 
boxes are located at the position corresponding to the studied locations. White lines rooted at 
those boxes indicate the magnetic field lines extrapolated from the magnetogram. 

to several weeks, the active region needs a certain mechanism to provide the plasma continuously. 
We note that the outflow region is localized at the edge of the active region, which means that a 
limited part of the active region is involved in the outflow. In contrast to this mass drain scenario, 
the extrapolated magnetic field lines rooted in the outflow region were connected to near the oppo¬ 
site edge of the active region according to the potential field calculation. The opposite side of the 
outflow region exhibit almost zero velocity, which indicates that the mass would accumulate from 
the outflow region. This leads to the picture that the outflow actually provides the active region 
with the plasma. However, the Doppler velocity map show a blueshifted pattern extending to the 
north west from the western outflow region, which may indicate that it is connected to far higher 
atmosphere. We must take into account the temporal evolution of the magnetic field in order to 
confirm the validity of these scenarios which will be studied in the near future. 
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6.3. Eastern and western outflow regions 


Here we discuss some implications for the coronal formation (i.e., heating) from the viewpoint 
of the outflows. The differences of derived quantities in those two outflow regions are listed in 
Table [TJ The topology of magnetic field lines can be inferred from the extrapolated field lines and 
the Doppler velocity map. We calculated the potential magnetic field from an MDI magnetogram 
taken during the EIS scan which started from 10:25:42UT, since its FOV is larger than that of the 
EIS scan used for the density diagnostics, and is large enough to include the entire active region. 
In order to confirm the connectivity of the magnetic field lines rooted at the studied locations, 
we drew projected field lines onto the intensity map of Fe xiv 264.78A as shown in Figure [l2j 
The outflow regions U1-U6 are indicated by white boxes. Note that since the intensity map was 
derived from the EIS scan which started from 10:25:42UT, we took into account the solar rotation 
to identify the locations of those boxes. The contours with orange (turquoise) indicate a magnetic 
field strength of +250 (—250) G in the MDI magnetogram. 


Two solid white lines trace coronal loops, therefore we regarded the topology of the eastern 
outflow region as closed, which can be also seen as a coherent pattern tracing the coronal loops 
in the Doppler velocity maps. Eour dashed white lines rooted at the western outflow region are 
connected to the opposite polarity around (x,y) = (—160", —150"), but the Doppler velocity maps 
clearly show that the blueshifted feature extends into the far west from which we suspected the 
topology of the western outflow region would be open. The closed loops rooted at the eastern 
outflow region are brighter than the open structures extending from the western outflow region by 
one order of magnitude. This might reflect the length of each structure in the sense that the upflow 
easily fills a closed loop while it flow s without obsta, c les in an open structure, which produces denser 
plasma in the closed loop. Note that ICulhane et al.l (j2fll4l ) suggested that the eastern outflow region 
is actually connected to the heliosphere through a two-step reconnection process. We may need 
further observations [e.g., statistical) in the future mission to clarify this point. 


As a consequence, it leads to the implication that the upflow from the bottom of the corona 
becomes dense in the closed loop because of the pressure balance between the corona and the 
transition region, which is consistent with our result that the electron density of EBW component 
was larger in the eastern region than in the western region (see Table [T]) . Although the difference 
in the electron density of the major component would not be insignificant, the relationship of the 
column depth (ie., larger /iMajor in the eastern outflow region than in the western outflow region) 
may represent that the eastern outflow region consists of more coronal loops than the western 
outflow region. 


We have evaluated mass leakage from the western out flow region in Section 16.21 A study on the 
first ionization potential (FIP) bias in the active region bv iBrooks fc WarrenI (j2fllll ) also suggested 
that the western outflow region connects to the slow solar wind. In contrast, the closed topology of 
the eastern outflow region may actually imply mass supply to the active region. If this is the case 
for a portion of the outflow region, it means that the outflow plays a crucial role in the coronal 
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Fig. 13.— Schematic picture of active region outflows. 

heating by supplying hot plasma into coronal loops. We suggest a possible picture in Figure [13] as 
a summary of this discussion. 


7. Summary 

The electron density of the outflow from the edges of NOAA AR10978 was measured by using 
an emission line pair Fe xiv 264.78A/274.20A. The upflow component was extracted from an 
enhanced blue wing (EBW) in Fe xiv line profiles through double-Gaussian fitting. We fitted those 
two Fe XIV emission lines simultaneously with a physical restriction that corresponding components 
in two emission lines must have the same Doppler velocity and thermal width, which previous EIS 
analysis on the density diagnostics have not been tried. The results were listed in Table [TJ 

The derived electron density for the major component (riMajor) and that for EBW component 
(nEBw) had opposite relationship in their magnitudes at the eastern and western outflow regions. 
There are several possibilities which cause the difference in the magnitude relationship between 
the east and west outflow region as follows. (1) The major component and EBW in Fe xiv line 
profiles are not directly related {e.g., superposition of structures along the line of sight). The 
electron density of EBW component just reflects the energy input amount. (2) The eastern outflow 
regions consist of the footpoints of corona loops extending to the north and connected to the 
opposite magnetic polarity around {x^y) = (—170",—70"), while longer coronal loops emanate in 
the western outflow regions and extend to the north west considering the appearance in Figure 0] 
The difference in length may influence the plasma density by the same driving mechanism for the 
outflow, since it is easier for the upflows in an open structure to flow without condensation than 
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for those in a closed loop. 

We also calculated the column depth for each component (/iMajor and Heew)- In the eastern 
region, h^BW was smaller than /iMajor by roughly one order of magnitude, which implies that the 
upflows possess only a small fraction (~ 0.1). Considering this implication with the result for the 
electron density (nEsw > IMajor)) h leads to a picture that the upflows may play a role in supplying 
hot plasma (logr[K] = 6.2-6.3) into coronal loops. On the other hand, in the western outflow 
region, the upflows have a larger volume by a factor of 5-6 than the plasma characterized by the 
major component, from which we consider the western outflow region as a structure composed of 
extending tubes with unidirectional upflows. 

We introduced a density diagnostics from a new point of view in Section [5l Electron density 
derived in our method is a function of Doppler velocity or wavelength (Equation 0]) , referred to as 
X-Ug diagram, which was found to be a good indicator of the electron density of minor components 
in a line prohle. The method has the advantage that it does not depend on any fitting model which 
might be ill-posed in some cases. Our aim was to evaluate the electron density of the outflow seen 
at the edge of the active region, and reinforce the result obtained in Section 01 

Using a density-sensitive emission line pair Ee xiv 264.78A/274.20A, we studied ne(A) by 
making X-rig diagrams at the active region core and the outflow regions. The increase in the 
diagrams was seen on the longer wavelength side for both structures, but we could not ascertain 
whether that behavior actually implies a physical situation at present. The diagrams for the 
active region core were flat around logUe [cm]“^ ~ 9.5, while those for the outflow regions exhibit 
some characteristic behaviors at the shorter wavelength side. They show a small hump around 
V = —110kms~^ in the eastern region (cuts 1 and 2 in Eigure[9]), and a decrease trend from 
logUe [cm“^] = 9.0 to logUe [cm“^] = 8.5 in a velocity scale of 100kms“^ in the western outflow 
region (cuts 3-5 in Eigure [9]) as seen in Eigure [TT] Thus we conhrmed the results obtained in 
Section 0] through our new method independent of the double-Gaussian fitting. 

As for the case where intermittent heating is responsible for the outflows, the duration of 
heating was crudely estimated to be longer than r = 500 s for the energy input of 10^^ erg (i.e., 
nanoflare) so that the density of upflows from the footpoints becomes compatible with that of the 
observed outflows. The electron density and column depth of the upflows in the eastern and western 
outflow regions were different, which was considered to be due to the magnetic structure above the 
outflow regions. Mass leakage occurs at the western outflow region (small ueew and large /iebw)- 
On the other hand, there is a possibility of the mass supply to active region loops at the eastern 
outflow region (large ueew and small /iebw); which may be related to the coronal heating process. 


Hinode is a Japanese mission developed and launched by ISAS/JAXA, with NAOJ as domestic 
partner and NASA and STEC (UK) as international partners. It is operated by these agencies in co¬ 
operation with ESA and NSC (Norway). We extend our gratitude to Hirohisa Hara and Toshifumi 
Shimizu who made a large number of insightful comments with their expertise in spectroscopy. We 
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Fig. 14.— Peak intensity of Fe xiv 264.78A vs. electron density deduced from the major component 
of Fe XIV 264.78A and 274.20A in the double Gaussian fitting. 


are also thankful to the GCOE program for proofreading/editing assistance. 
Facilities: Hinode. 


A. Fe xiv 264.78A intensity and electron density 

We can see a clear positive correlation between peak intensity and electron density in the 
intensity range larger than /Major = 2.0 x 10^ ergcm“^ s“^ sr“^ A ^ in Figure O (indicated by a 
vertical dashed line) while the plot is more scattered below that intensity. Not only the photon 
noise contributes to this large degree of uncertainty, but also unidentified blended emission lines 
could do so. Therefore we analyzed the data points with /Major larger than the value which the 
vertical dashed line indicates. 


B. Uncertainty in Si vii density 

Since the electron density is not the same for emission lines with different formation temper¬ 
ature, there is an uncertainty in nsivii which cannot be determined from the data used in this 
analysis. In order to evaluate the error in the electron density derived for Fe xiv (npexiv) coming 
from this uncertainty, we remove the blending Si vii at Fe xiv 274.20A in three cases for resivin 10®, 
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Fig. 15.— (a) Scatter plots of derived Fe xiv density (npexiv) for different electron density of 
Si VII (nsivii). Horizontal axis indicates npexiv of EBW component derived by assuming nsivu = 
10® cm~^. Vertical axis indicates npexiv of EBW component derived by assuming nsivo = 10® cm~® 
(blue) and 10^® cm“® (red), (b) The same data as in panel (a) but the horizontal axis indicates 
a relative difference Anpexiv/fiPexiv, where Anpexiv is a difference of npexiv for different nsivu (10® 
and 10^®cm~^) measured from the case for nsivu = 10®cm“®. 

10®, and 10^®cm“®, and derived npexiv for each case. Panel (a) in Figure fT5] shows scatter plots 
for the electron density of EBW component within the entire western outflow region derived for 
the case nsivu = 10® cm“® vs. 10® cm“® (10^®cm“®) in blue (red). The npexiv of EBW component 
derived by assuming nsivu = 10® (10^®)cm“® becomes smaller (larger). Panel (b) in Figure [15] 
shows those relative differences Anpexiv/ii-FexiV) where Anpexiv is a difference of npexiv for different 
rrsivii (10® and 10^®cm“^) measured from the case for nsivu = 10®cm“®. Colors {red and blue) are 
the same as in panel (a). Solid and Dashed histograms indicate the western outflow region (the 
white dashed box in Figure [5|) and for the entire field of view, respectively. These relative differ¬ 
ences were calculated in log scale. The histograms show that the error coming from the difference 
of nsivii does not exceed 5%. It means that the error is around at most for the density 

range 10® cm“® < ne < 10^® cm“®, and roughly becomes a factor of 3 (f.e., comparable to the error 
originated in the photon noise). 
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